Direct impact of COVID-19 vaccination in Chile: averted cases, hospitalizations, ICU admissions, and deaths.

BACKGROUND: Chile rapidly implemented an extensive COVID-19 vaccination campaign, deploying a diversity of vaccines with a strategy that prioritized the elderly and individuals with comorbidities. This study aims to assess the direct impact of vaccination on the number of COVID-19 related cases, hospital admissions, ICU admissions and deaths averted during the first year and a half of the campaign. METHODS: Via Chile's transparency law, we obtained access to weekly event counts categorized by vaccination status and age. Integrating this data with publicly available census and vaccination coverage information, we conducted a comparative analysis of weekly incidence rates between vaccinated and unvaccinated groups from December 20, 2020 to July 2, 2022 to estimate the direct impact of vaccination in terms of the number of cases, hospitalizations, ICU admissions and deaths averted, using an approach that avoids the need to explicitly specify the effectiveness of each vaccine deployed. RESULTS: We estimated that, from December 20, 2020 to July 2, 2022 the vaccination campaign directly prevented 1,030,648 (95% Confidence Interval: 1,016,975-1,044,321) cases, 268,784 (95% CI: 264,524-273,045) hospitalizations, 85,830 (95% CI: 83,466-88,194) ICU admissions and 75,968 (95% CI: 73,909-78,028) deaths related to COVID-19 among individuals aged 16 years and older. This corresponds to a reduction of 26% of cases, 66% of hospital admissions, 70% of ICU admissions and 67% of deaths compared to a scenario without vaccination. Individuals 55 years old or older represented 67% of hospitalizations, 73% of ICU admissions and 89% of deaths related to COVID-19 prevented. CONCLUSIONS: This study highlights the role of Chile's vaccination campaign in reducing COVID-19 disease burden, with the most substantial reductions observed in severe outcomes.

Unveiling abundance-dependent metabolic phenotypes of microbial communities.

IMPORTANCE: In nature, organisms live in communities and not as isolated species, and their interactions provide a source of resilience to environmental disturbances. Despite their importance in ecology, human health, and industry, understanding how organisms interact in different environments remains an open question. In this work, we provide a novel approach that, only using genomic information, studies the metabolic phenotype exhibited by communities, where the exploration of suboptimal growth flux distributions and the composition of a community allows to unveil its capacity to respond to environmental changes, shedding light of the degrees of metabolic plasticity inherent to the community.

DNA sequencing in the classroom: complete genome sequence of two earwig (Dermaptera; Insecta) species.

BACKGROUND: Despite representing the largest fraction of animal life, the number of insect species whose genome has been sequenced is barely in the hundreds. The order Dermaptera (the earwigs) suffers from a lack of genomic information despite its unique position as one of the basally derived insect groups and its importance in agroecosystems. As part of a national educational and outreach program in genomics, a plan was formulated to engage the participation of high school students in a genome sequencing project. Students from twelve schools across Chile were instructed to capture earwig specimens in their geographical area, to identify them and to provide material for genome sequencing to be carried out by themselves in their schools. RESULTS: The school students collected specimens from two cosmopolitan earwig species: Euborellia annulipes (Fam. Anisolabididae) and Forficula auricularia (Fam. Forficulidae). Genomic DNA was extracted and, with the help of scientific teams that traveled to the schools, was sequenced using nanopore sequencers. The sequence data obtained for both species was assembled and annotated. We obtained genome sizes of 1.18 Gb (F. auricularia) and 0.94 Gb (E. annulipes) with the number of predicted protein coding genes being 31,800 and 40,000, respectively. Our analysis showed that we were able to capture a high percentage (≥ 93%) of conserved proteins indicating genomes that are useful for comparative and functional analysis. We were also able to characterize structural elements such as repetitive sequences and non-coding RNA genes. Finally, functional categories of genes that are overrepresented in each species suggest important differences in the process underlying the formation of germ cells, and modes of reproduction between them, features that are one of the distinguishing biological properties that characterize these two distant families of Dermaptera. CONCLUSIONS: This work represents an unprecedented instance where the scientific and lay community have come together to collaborate in a genome sequencing project. The versatility and accessibility of nanopore sequencers was key to the success of the initiative. We were able to obtain full genome sequences of two important and widely distributed species of insects which had not been analyzed at this level previously. The data made available by the project should illuminate future studies on the Dermaptera.

Plant and buffer state estimation for networked predictive control over multiple erasure channels.

This paper investigates optimal estimation for networked control systems composed of linear, time-variant, and multi-input multi-output (MIMO) plants that communicate with the controller over multiple independent erasure channels. In the controller-to-plant link, we cope with data losses by considering a networked predictive control strategy in which the controller not only sends the current control inputs to the actuators, but also sequences of predicted control inputs at every time step. These are stored in buffers located at the corresponding actuator side, with the aim of being used in case of data loss. Additionally, for losses occurring in the plant-to-controller link, we assume that the controller contains a data-dropout compensation mechanism capable of replacing the missing measurements with their optimal estimates. It is assumed that acknowledgment packages are sent from the plant to the controller, to indicate if the previous sequence was received or not. For this setup, we compute optimal time-varying estimates of the plant and buffer state, and for the lost measurements as well. Finally, we provide the optimal stationary estimator, which serves as a simpler alternative to the time-varying one.

DNA sequencing in the classroom: complete genome sequence of two earwig (Dermaptera; Insecta) species

BACKGROUND: Despite representing the largest fraction of animal life, the number of insect species whose genome has been sequenced is barely in the hundreds. The order Dermaptera (the earwigs) suffers from a lack of genomic information despite its unique position as one of the basally derived insect groups and its importance in agroecosystems. As part of a national educational and outreach program in genomics, a plan was formulated to engage the participation of high school students in a genome sequencing project. Students from twelve schools across Chile were instructed to capture earwig specimens in their geographical area, to identify them and to provide material for genome sequencing to be carried out by themselves in their schools. RESULTS: The school students collected specimens from two cosmopolitan earwig species: Euborellia annulipes (Fam. Anisolabididae) and Forficula auricularia (Fam. Forficulidae). Genomic DNA was extracted and, with the help of scientific teams that traveled to the schools, was sequenced using nanopore sequencers. The sequence data obtained for both species was assembled and annotated. We obtained genome sizes of 1.18 Gb (F. auricularia) and 0.94 Gb (E. annulipes) with the number of predicted protein coding genes being 31,800 and 40,000, respectively. Our analysis showed that we were able to capture a high percentage (≥ 93%) of conserved proteins indicating genomes that are useful for comparative and functional analysis. We were also able to characterize structural elements such as repetitive sequences and non-coding RNA genes. Finally, functional categories of genes that are overrepresented in each species suggest important differences in the process underlying the formation of germ cells, and modes of reproduction between them, features that are one of the distinguishing biological properties that characterize these two distant families of Dermaptera. CONCLUSIONS: This work represents an unprecedented instance where the scientific and lay community have come together to collaborate in a genome sequencing project. The versatility and accessibility of nanopore sequencers was key to the success of the initiative. We were able to obtain full genome sequences of two important and widely distributed species of insects which had not been analyzed at this level previously. The data made available by the project should illuminate future studies on the Dermaptera.

Genome sequencing and transcriptomic analysis of the Andean killifish Orestias ascotanensis reveals adaptation to high-altitude aquatic life.

Orestias ascotanensis (Cyprinodontidae) is a teleost pupfish endemic to springs feeding into the Ascotan saltpan in the Chilean Altiplano (3,700 m.a.s.l.) and represents an opportunity to study adaptations to high-altitude aquatic environments. We have de novo assembled the genome of O. ascotanensis at high coverage. Comparative analysis of the O. ascotanensis genome showed an overall process of contraction, including loss of genes related to G-protein signaling, chemotaxis and signal transduction, while there was expansion of gene families associated with microtubule-based movement and protein ubiquitination. We identified 818 genes under positive selection, many of which are involved in DNA repair. Additionally, we identified novel and conserved microRNAs expressed in O. ascotanensis and its closely-related species, Orestias gloriae. Our analysis suggests that positive selection and expansion of genes that preserve genome stability are a potential adaptive mechanism to cope with the increased solar UV radiation to which high-altitude animals are exposed to.

Plant ecological genomics at the limits of life in the Atacama Desert.

The Atacama Desert in Chile-hyperarid and with high-ultraviolet irradiance levels-is one of the harshest environments on Earth. Yet, dozens of species grow there, including Atacama-endemic plants. Herein, we establish the Talabre-Lejía transect (TLT) in the Atacama as an unparalleled natural laboratory to study plant adaptation to extreme environmental conditions. We characterized climate, soil, plant, and soil-microbe diversity at 22 sites (every 100 m of altitude) along the TLT over a 10-y period. We quantified drought, nutrient deficiencies, large diurnal temperature oscillations, and pH gradients that define three distinct vegetational belts along the altitudinal cline. We deep-sequenced transcriptomes of 32 dominant plant species spanning the major plant clades, and assessed soil microbes by metabarcoding sequencing. The top-expressed genes in the 32 Atacama species are enriched in stress responses, metabolism, and energy production. Moreover, their root-associated soils are enriched in growth-promoting bacteria, including nitrogen fixers. To identify genes associated with plant adaptation to harsh environments, we compared 32 Atacama species with the 32 closest sequenced species, comprising 70 taxa and 1,686,950 proteins. To perform phylogenomic reconstruction, we concatenated 15,972 ortholog groups into a supermatrix of 8,599,764 amino acids. Using two codon-based methods, we identified 265 candidate positively selected genes (PSGs) in the Atacama plants, 64% of which are located in Pfam domains, supporting their functional relevance. For 59/184 PSGs with an Arabidopsis ortholog, we uncovered functional evidence linking them to plant resilience. As some Atacama plants are closely related to staple crops, these candidate PSGs are a "genetic goldmine" to engineer crop resilience to face climate change.

Nutrient Scarcity in a New Defined Medium Reveals Metabolic Resistance to Antibiotics in the Fish Pathogen Piscirickettsia salmonis.

Extensive use of antibiotics has been the primary treatment for the Salmonid Rickettsial Septicemia, a salmonid disease caused by the bacterium Piscirickettsia salmonis. Occurrence of antibiotic resistance has been explored in various P. salmonis isolates using different assays; however, P. salmonis is a nutritionally demanding intracellular facultative pathogen; thus, assessing its antibiotic susceptibility with standardized and validated protocols is essential. In this work, we studied the pathogen response to antibiotics using a genomic, a transcriptomic, and a phenotypic approach. A new defined medium (CMMAB) was developed based on a metabolic model of P. salmonis. CMMAB was formulated to increase bacterial growth in nutrient-limited conditions and to be suitable for performing antibiotic susceptibility tests. Antibiotic resistance was evaluated based on a comprehensive search of antibiotic resistance genes (ARGs) from P. salmonis genomes. Minimum inhibitory concentration assays were conducted to test the pathogen susceptibility to antibiotics from drug categories with predicted ARGs. In all tested P. salmonis strains, resistance to erythromycin, ampicillin, penicillin G, streptomycin, spectinomycin, polymyxin B, ceftazidime, and trimethoprim was medium-dependent, showing resistance to higher antibiotic concentrations in the CMMAB medium. The mechanism for antibiotic resistance to ampicillin in the defined medium was further explored and was proven to be associated to a decrease in the bacterial central metabolism, including the TCA cycle, the pentose-phosphate pathway, energy production, and nucleotide metabolism, and it was not associated with decreased growth rate of the bacterium or with the expression of any predicted ARG. Our results suggest that nutrient scarcity plays a role in the bacterial antibiotic resistance, protecting against the detrimental effects of antibiotics, and thus, we propose that P. salmonis exhibits a metabolic resistance to ampicillin when growing in a nutrient-limited medium.

Global Proteomic Profiling of Piscirickettsia salmonis and Salmon Macrophage-Like Cells during Intracellular Infection.

Piscirickettsiasalmonis is an intracellular bacterial fish pathogen that causes piscirickettsiosis, a disease with numerous negative impacts in the Chilean salmon farming industry. Although transcriptomic studies of P. salmonis and its host have been performed, dual host-pathogen proteomic approaches during infection are still missing. Considering that gene expression does not always correspond with observed phenotype, and bacteriological culture studies inadequately reflect infection conditions, to improve the existing knowledge for the pathogenicity of P. salmonis, we present here a global proteomic profiling of Salmon salar macrophage-like cell cultures infected with P. salmonis LF-89. The proteomic analyses identified several P. salmonis proteins from two temporally different stages of macrophages infection, some of them related to key functions for bacterial survival in other intracellular pathogens. Metabolic differences were observed in early-stage infection bacteria, compared to late-stage infections. Virulence factors related to membrane, lipopolysaccharide (LPS) and surface component modifications, cell motility, toxins, and secretion systems also varied between the infection stages. Pilus proteins, beta-hemolysin, and the type VI secretion system (T6SS) were characteristic of the early-infection stage, while fimbria, upregulation of 10 toxins or effector proteins, and the Dot/Icm type IV secretion system (T4SS) were representative of the late-infection stage bacteria. Previously described virulence-related genes in P. salmonis plasmids were identified by proteomic assays during infection in SHK-1 cells, accompanied by an increase of mobile-related elements. By comparing the infected and un-infected proteome of SHK-1 cells, we observed changes in cellular and redox homeostasis; innate immune response; microtubules and actin cytoskeleton organization and dynamics; alteration in phagosome components, iron transport, and metabolism; and amino acids, nucleoside, and nucleotide metabolism, together with an overall energy and ATP production alteration. Our global proteomic profiling and the current knowledge of the P. salmonis infection process allowed us to propose a model of the macrophage-P. salmonis interaction.

Identification of SNPs and InDels associated with berry size in table grapes integrating genetic and transcriptomic approaches.

BACKGROUND: Berry size is considered as one of the main selection criteria in table grapes breeding programs, due to the consumer preferences. However, berry size is a complex quantitive trait under polygenic control, and its genetic determination of berry weight is not yet fully understood. The aim of this work was to perform marker discovery using a transcriptomic approach, in order to identify and characterize SNP and InDel markers associated with berry size in table grapes. We used an integrative analysis based on RNA-Seq, SNP/InDel search and validation on table grape segregants and varieties with different genetic backgrounds. RESULTS: Thirty SNPs and eight InDels were identified using a transcriptomic approach (RNA-Seq). These markers were selected from SNP/InDel found among segregants from a Ruby x Sultanina population with contrasting phenotypes for berry size. The set of 38 SNP and InDel markers was distributed in eight chromosomes. Genotype-phenotype association analyses were performed using a set of 13 RxS segregants and 41 table grapes varieties with different genetic backgrounds during three seasons. The results showed several degrees of association of these markers with berry size (10.2 to 30.7%) as other berry-related traits such as length and width. The co-localization of SNP and /or InDel markers and previously reported QTLs and candidate genes associated with berry size were analysed. CONCLUSIONS: We identified a set of informative and transferable SNP and InDel markers associated with berry size. Our results suggest the suitability of SNPs and InDels as candidate markers for berry weight in seedless table grape breeding. The identification of genomic regions associated with berry weight in chromosomes 8, 15 and 17 was achieved with supporting evidence derived from a transcriptome experiment focused on SNP/InDel search, as well as from a QTL-linkage mapping approach. New regions possibly associated with berry weight in chromosomes 3, 6, 9 and 14 were identified.

Development of a small panel of SNPs to infer ancestry in Chileans that distinguishes Aymara and Mapuche components.

BACKGROUND: Current South American populations trace their origins mainly to three continental ancestries, i.e. European, Amerindian and African. Individual variation in relative proportions of each of these ancestries may be confounded with socio-economic factors due to population stratification. Therefore, ancestry is a potential confounder variable that should be considered in epidemiologic studies and in public health plans. However, there are few studies that have assessed the ancestry of the current admixed Chilean population. This is partly due to the high cost of genome-scale technologies commonly used to estimate ancestry. In this study we have designed a small panel of SNPs to accurately assess ancestry in the largest sampling to date of the Chilean mestizo population (n = 3349) from eight cities. Our panel is also able to distinguish between the two main Amerindian components of Chileans: Aymara from the north and Mapuche from the south. RESULTS: A panel of 150 ancestry-informative markers (AIMs) of SNP type was selected to maximize ancestry informativeness and genome coverage. Of these, 147 were successfully genotyped by KASPar assays in 2843 samples, with an average missing rate of 0.012, and a 0.95 concordance with microarray data. The ancestries estimated with the panel of AIMs had relative high correlations (0.88 for European, 0.91 for Amerindian, 0.70 for Aymara, and 0.68 for Mapuche components) with those obtained with AXIOM LAT1 array. The country's average ancestry was 0.53 ± 0.14 European, 0.04 ± 0.04 African, and 0.42 ± 0.14 Amerindian, disaggregated into 0.18 ± 0.15 Aymara and 0.25 ± 0.13 Mapuche. However, Mapuche ancestry was highest in the south (40.03%) and Aymara in the north (35.61%) as expected from the historical location of these ethnic groups. We make our results available through an online app and demonstrate how it can be used to adjust for ancestry when testing association between incidence of a disease and nongenetic risk factors. CONCLUSIONS: We have conducted the most extensive sampling, across many different cities, of current Chilean population. Ancestry varied significantly by latitude and human development. The panel of AIMs is available to the community for estimating ancestry at low cost in Chileans and other populations with similar ancestry.

Genome-scale metabolic models of Microbacterium species isolated from a high altitude desert environment.

The Atacama Desert is the most arid desert on Earth, focus of important research activities related to microbial biodiversity studies. In this context, metabolic characterization of arid soil bacteria is crucial to understand their survival strategies under extreme environmental stress. We investigated whether strain-specific features of two Microbacterium species were involved in the metabolic ability to tolerate/adapt to local variations within an extreme desert environment. Using an integrative systems biology approach we have carried out construction and comparison of genome-scale metabolic models (GEMs) of two Microbacterium sp., CGR1 and CGR2, previously isolated from physicochemically contrasting soil sites in the Atacama Desert. Despite CGR1 and CGR2 belong to different phylogenetic clades, metabolic pathways and attributes are highly conserved in both strains. However, comparison of the GEMs showed significant differences in the connectivity of specific metabolites related to pH tolerance and CO2 production. The latter is most likely required to handle acidic stress through decarboxylation reactions. We observed greater GEM connectivity within Microbacterium sp. CGR1 compared to CGR2, which is correlated with the capacity of CGR1 to tolerate a wider pH tolerance range. Both metabolic models predict the synthesis of pigment metabolites (ß-carotene), observation validated by HPLC experiments. Our study provides a valuable resource to further investigate global metabolic adaptations of bacterial species to grow in soils with different abiotic factors within an extreme environment.

High-Throughput Single Nucleotide Polymorphism (SNP) Discovery and Validation Through Whole-Genome Resequencing in Nile Tilapia (Oreochromis niloticus).

Nile tilapia (Oreochromis niloticus) is the second most important farmed fish in the world and a sustainable source of protein for human consumption. Several genetic improvement programs have been established for this species in the world. Currently, the estimation of genetic merit of breeders is typically based on genealogical and phenotypic information. Genome-wide information can be exploited to efficiently incorporate traits that are difficult to measure into the breeding goal. Thus, single nucleotide polymorphisms (SNPs) are required to investigate phenotype-genotype associations and determine the genomic basis of economically important traits. We performed de novo SNP discovery in three different populations of farmed Nile tilapia. A total of 29.9 million non-redundant SNPs were identified through Illumina (HiSeq 2500) whole-genome resequencing of 326 individual samples. After applying several filtering steps, including removing SNP based on genotype and site quality, presence of Mendelian errors, and non-unique position in the genome, a total of 50,000 high-quality SNPs were selected for the development of a custom Illumina BeadChip SNP panel. These SNPs were highly informative in the three populations analyzed showing between 43,869 (94%) and 46,139 (99%) SNPs in Hardy-Weinberg Equilibrium; 37,843 (76%) and 45,171(90%) SNPs with a minor allele frequency (MAF) higher than 0.05; and 43,450 (87%) and 46,570 (93%) SNPs with a MAF higher than 0.01. The 50K SNP panel developed in the current work will be useful for the dissection of economically relevant traits, enhancing breeding programs through genomic selection, as well as supporting genetic studies in farmed populations of Nile tilapia using dense genome-wide information.

Development of a small panel of SNPs to infer ancestry in Chileans that distinguishes Aymara and Mapuche components

BACKGROUND: Current South American populations trace their origins mainly to three continental ancestries, i.e. European, Amerindian and African. Individual variation in relative proportions of each of these ancestries may be confounded with socio-economic factors due to population stratification. Therefore, ancestry is a potential confounder variable that should be considered in epidemiologic studies and in public health plans. However, there are few studies that have assessed the ancestry of the current admixed Chilean population. This is partly due to the high cost of genome-scale technologies commonly used to estimate ancestry. In this study we have designed a small panel of SNPs to accurately assess ancestry in the largest sampling to date of the Chilean mestizo population (n = 3349) from eight cities. Our panel is also able to distinguish between the two main Amerindian components of Chileans: Aymara from the north and Mapuche from the south. RESULTS: A panel of 150 ancestry-informative markers (AIMs) of SNP type was selected to maximize ancestry informativeness and genome coverage. Of these, 147 were successfully genotyped by KASPar assays in 2843 samples, with an average missing rate of 0.012, and a 0.95 concordance with microarray data. The ancestries estimated with the panel of AIMs had relative high correlations (0.88 for European, 0.91 for Amerindian, 0.70 for Aymara, and 0.68 for Mapuche components) with those obtained with AXIOM LAT1 array. The country's average ancestry was 0.53 ± 0.14 European, 0.04 ± 0.04 African, and 0.42 ± 0.14 Amerindian, disaggregated into 0.18 ± 0.15 Aymara and 0.25 ± 0.13 Mapuche. However, Mapuche ancestry was highest in the south (40.03%) and Aymara in the north (35.61%) as expected from the historical location of these ethnic groups. We make our results available through an online app and demonstrate how it can be used to adjust for ancestry when testing association between incidence of a disease and nongenetic risk factors. CONCLUSIONS: We have conducted the most extensive sampling, across many different cities, of current Chilean population. Ancestry varied significantly by latitude and human development. The panel of AIMs is available to the community for estimating ancestry at low cost in Chileans and other populations with similar ancestry.

Piscirickettsia salmonis Cryptic Plasmids: Source of Mobile DNA and Virulence Factors.

Four large cryptic plasmids were identified in the salmon pathogen Piscirickettsia salmonis reference strain LF-89. These plasmids appeared highly novel, with less than 7% nucleotidic identity to the nr plasmid database. Plasmid copy number analysis revealed that they are harbored in chromosome equivalent ratios. In addition to plasmid-related genes (plasmidial autonomous replication, partitioning, maintenance, and mobilization genes), mobile genetic elements such as transposases, integrases, and prophage sequences were also identified in P. salmonis plasmids. However, bacterial lysis was not observed upon the induction of prophages. A total of twelve putative virulence factors (VFs) were identified, in addition to two global transcriptional regulators, the widely conserved CsrA protein and the regulator Crp/Fnr. Eleven of the putative VFs were overexpressed during infection in two salmon-derived cellular infection models, supporting their role as VFs. The ubiquity of these plasmids was also confirmed by sequence similarity in the genomes of other P. salmonis strains. The ontology of P. salmonis plasmids suggests a role in bacterial fitness and adaptation to the environment as they encode proteins related to mobilization, nutrient transport and utilization, and bacterial virulence. Further functional characterization of P. salmonis plasmids may improve our knowledge regarding virulence and mobile elements in this intracellular pathogen.

Fine Mapping Using Whole-Genome Sequencing Confirms Anti-Müllerian Hormone as a Major Gene for Sex Determination in Farmed Nile Tilapia (Oreochromis niloticus L.).

Nile tilapia (Oreochromis niloticus) is one of the most cultivated and economically important species in world aquaculture. Intensive production promotes the use of monosex animals, due to an important dimorphism that favors male growth. Currently, the main mechanism to obtain all-male populations is the use of hormones in feeding during larval and fry phases. Identifying genomic regions associated with sex determination in Nile tilapia is a research topic of great interest. The objective of this study was to identify genomic variants associated with sex determination in three commercial populations of Nile tilapia. Whole-genome sequencing of 326 individuals was performed, and a total of 2.4 million high-quality bi-allelic single nucleotide polymorphisms (SNPs) were identified after quality control. A genome-wide association study (GWAS) was conducted to identify markers associated with the binary sex trait (males = 1; females = 0). A mixed logistic regression GWAS model was fitted and a genome-wide significant signal comprising 36 SNPs, spanning a genomic region of 536 kb in chromosome 23 was identified. Ten out of these 36 genetic variants intercept the anti-Müllerian (Amh) hormone gene. Other significant SNPs were located in the neighboring Amh gene region. This gene has been strongly associated with sex determination in several vertebrate species, playing an essential role in the differentiation of male and female reproductive tissue in early stages of development. This finding provides useful information to better understand the genetic mechanisms underlying sex determination in Nile tilapia.

Expanding Metabolic Capabilities Using Novel Pathway Designs: Computational Tools and Case Studies.

Design and selection of efficient metabolic pathways is critical for the success of metabolic engineering endeavors. Convenient pathways should not only produce the target metabolite in high yields but also are required to be thermodynamically feasible under production conditions, and to prefer efficient enzymes. To support the design and selection of such pathways, different computational approaches have been proposed for exploring the feasible pathway space under many of the above constraints. In this review, an overview of recent constraint-based optimization frameworks for metabolic pathway prediction, as well as relevant pathway engineering case studies that highlight the importance of rational metabolic designs is presented. Despite the availability and suitability of in silico design tools for metabolic pathway engineering, scarce-although increasing-application of computational outcomes is found. Finally, challenges and limitations hindering the broad adoption and successful application of these tools in metabolic engineering projects are discussed.

Creation and analysis of biochemical constraint-based models using the COBRA Toolbox v.3.0.

Constraint-based reconstruction and analysis (COBRA) provides a molecular mechanistic framework for integrative analysis of experimental molecular systems biology data and quantitative prediction of physicochemically and biochemically feasible phenotypic states. The COBRA Toolbox is a comprehensive desktop software suite of interoperable COBRA methods. It has found widespread application in biology, biomedicine, and biotechnology because its functions can be flexibly combined to implement tailored COBRA protocols for any biochemical network. This protocol is an update to the COBRA Toolbox v.1.0 and v.2.0. Version 3.0 includes new methods for quality-controlled reconstruction, modeling, topological analysis, strain and experimental design, and network visualization, as well as network integration of chemoinformatic, metabolomic, transcriptomic, proteomic, and thermochemical data. New multi-lingual code integration also enables an expansion in COBRA application scope via high-precision, high-performance, and nonlinear numerical optimization solvers for multi-scale, multi-cellular, and reaction kinetic modeling, respectively. This protocol provides an overview of all these new features and can be adapted to generate and analyze constraint-based models in a wide variety of scenarios. The COBRA Toolbox v.3.0 provides an unparalleled depth of COBRA methods.

Whole Genome Sequence, Variant Discovery and Annotation in Mapuche-Huilliche Native South Americans.

Whole human genome sequencing initiatives help us understand population history and the basis of genetic diseases. Current data mostly focuses on Old World populations, and the information of the genomic structure of Native Americans, especially those from the Southern Cone is scant. Here we present annotation and variant discovery from high-quality complete genome sequences of a cohort of 11 Mapuche-Huilliche individuals (HUI) from Southern Chile. We found approximately 3.1 × 106 single nucleotide variants (SNVs) per individual and identified 403,383 (6.9%) of novel SNVs events. Analyses of large-scale genomic events detected 680 copy number variants (CNVs) and 4,514 structural variants (SVs), including 398 and 1,910 novel events, respectively. Global ancestry composition of HUI genomes revealed that the cohort represents a sample from a marginally admixed population from the Southern Cone, whose main genetic component derives from Native American ancestors. Additionally, we found that HUI genomes contain variants in genes associated with 5 of the 6 leading causes of noncommunicable diseases in Chile, which may have an impact on the risk of prevalent diseases in Chilean and Amerindian populations. Our data represents a useful resource that can contribute to population-based studies and for the design of early diagnostics or prevention tools for Native and admixed Latin American populations.

Bacterial communities associated to Chilean altiplanic native plants from the Andean grasslands soils.

The rhizosphere is considered the primary place for soil microbiome differentiation and plays a key role in plant survival, especially for those subjected to environmental stress. Using high-throughput sequencing of the 16S rRNA gene, we analyzed and compared soil bacterial communities associated to four of the most abundant high altitude native plant species of the Chilean Andean grasslands. We examined three soil compartments: the rhizosphere (bacteria firmly attached to the roots), the rhizosphere-surrounding soil (bacteria loosely attached to the roots) and the bulk soil (plant-free soil). The rhizosphere microbiome was in all cases the least diverse, exposing that the bulk soil was a more complex environment. Taxonomic analysis revealed an abrupt change between the rhizosphere and the rest of the non-rhizospheric soils. Thus, while rhizobacterial communities were enriched in Proteobacteria (mainly Alphaproteobacteria), Actinobacteria (mostly Blastocatellia) dominated in bulk soils. Finally, we detected certain taxonomic rhizosphere signatures, which could be attributed to a particular genotype. Overall, our results indicate that the thin layer of soil surrounding the roots constitute a distinctive soil environment. This study contributes to expand the knowledge about soil bacterial communities in the Chilean highlands and takes the first step to understand the processes that might lead to the rhizosphere differentiation in that area.